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Active layer model for wall-bounded turbulence 

By M. T. Landahl 1 J. Kim 2 and P. R. Spalart 2 


The active-layer model for wall-bounded turbulence hypothesizes that the non- 
linear terms are large only in a thin layer near the wall, and hence the turbulence 
in the region outside the active inner layer can be modeled as a linear fluctuating 
flow driven by the active layer. This hypothesis is tested using data obtained from 
a direct simulation of turbulent channel flow. It is found that the nonlinear effects 
are the strongest near the wall with a maximum at around y + = 20 and, outside 
the near- wall region, these involve primarily the cascading mechanism leading to 
dissipation. 


1. Introduction 

Laboratory experiments and numerical simulation have shown that the turbulent 
activity in wall-bounded turbulence is the highest in the immediate neighborhood 
of the wall, in the viscous and buffer layers. Therefore, nonlinear effects may be 
expected to be strongest in this region. A possible model for the turbulent field 
might therefore be to consider the turbulence in the region outside an active inner 
layer as a linear fluctuating flow driven by the active layer (Fig. 1). 

2. Analysis 

We subdivide the flow field in a parallel mean flow, U(y)6n , and a fluctuating 

field, u,v,w,p. By elimination of the pressure from the momentum equations, 

making use of the continuity equations in the process, the following equation for 
the v- component is found (Landahl, 1967): 

Los(v) — q (1) 

where Los is the (space-time) Orr-Sommerfeld operator, 

Los(v) = ( d/dt + Ud/dx)V 2 v - U"dv/dx - uV 4 v (2) 

and q contains all the nonlinear terms, 

q = V 2 T 2 - d 2 Ti/dxidx 2 (3) 

where 

Ti = drij/dxj (4) 
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Figure 1. A model of wall-bounded turbulence consist of active non-linear inner 
layer and a linear fluctuating outer layer driven by the active layer. 


with 

Tij = -UiUj (5) 

The solution of (1) assuming q to be given may be determined by applying Fourier 
transformation in x , z and f, which yields 

(17 - c)(v" - k 2 v) - U"v - ( l/ik x ){d 2 /dy 2 - k 2 ) 2 v = q/ik x (6) 

where the caret denotes a Fourier-transformed quantity, k x and k z , denote the x- 
and r- components of the wave number k, respectively, and c = u>/k x . 

3. Calculation of the Reynolds stress 

The basic hypotheses in the active-layer model are that the nonlinear terms are 
large only in a thin region near the wall and that the flow in the region outside 
this layer may be found by considering q as a given and driving an outer linear 
fluctuating flow field. From the statistics of q one may then determine any desired 
statistical quantity such as the mean Reynolds shear stress, 

r = —p < uv > (7) 

By applying Fourier transform to the component equations and using Parseval’s 
formula relating quadratic mean quantities to their transforms, one finds, upon 
neglecting viscous terms, that the wave number spectrum, S T , for the Reynolds 
shear stress, is related to the spectrum for the v-fluctuations through 

s ;=( fc 2 / fc ) 2 ( d/dy)(U'S c v ) + (k x /k ) 2 U" SI ( 8 ) 

where 5' denotes the spectrum of v in a reference frame convected with the mean 
velocity I7(y), and the prime denotes differentiation with respect to y. From (6), it 
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follows that for small viscosity the v- fluctuations may receive their largest contri- 
butions from the region near k x = 0, in which case the second term in (8) will be 
negligible, and with k = k z one finds after integration that 

Sr = U'S C V (9) 

a result that may be shown to be consistent with Prandtl’s mixing-length hypothesis. 

4. Application to channel flow turbulence 

The basic hypothesis on which the present model is based, namely, that the non- 
linearities are important only in a thin region near the wall, may be tested with 
the aid of the numerical turbulence simulations. A fairly extensive data base for 
turbulent channel flow has been generated in the NASA-Ames numerical simula- 
tions, although limited to fairly low Reynolds numbers. From the computed data, 
the nonlinear terms of (1) may be extracted and appropriate statistical quantities 
determined. Since the detailed temporal evolution is not so easily accessible because 
of the way the data are stored, for a preliminary assessment of the soundness of 
the basic hypothesis of the present model, the spatial mean and power spectrum 
of a single time realization was determined from the computed velocity field. The 
results are shown in Figs. 2 and 3. In Fig. 2, the root-mean- square value of q is 
presented (arbitrary scale). It has a maximum at around y/h = 0.1 ( h — channel 
half width), which corresponds to y+ = 18, and drops off to a value of about one 
tenth of the maximum towards the center of the channel. The spectra shown in 
Fig. 3 for different distances from the wall are found to be fairly flat in the wall 
and buffer regions, with the cutoff in k z at a higher value than for k x , reflecting 
the dominance of high streamwise elongation of the dominating structures. The 
spectrum at y + = 76 (Figs. 3g and 3h), however, which is well outside the buffer 
layer, has its main contribution from the high wave- number range (both in k x and 
k z ). This probably reflects the nonlinear cascading mechanism involved in dissipa- 
tion. The production of Reynolds stress resulting from this spectrum is not large, 
however, since it receives its main contributions from the low A;*, -end. The part of 
the ^-spectrum responsible for the production can therefore be expected to be even 
more concentrated in the near-wall region. 

5. Conclusion 

The model proposed for wall-bounded turbulence, namely, that the nonlinear 
driving of the turbulent fluctuations is concentrated in the near-wall region, in 
the viscous and buffer regions, has been tentatively examined using data from the 
NASA-Ames channel-flow simulations. The preliminary findings are that the non- 
linear effects are the strongest near the wall, with a maximum at around y + = 20, 
and that outside the near-wall region they involve primarily the cascading mecha- 
nism leading to dissipation. The mean and the spectra were obtained from a single 
time realization; for a more general treatment, one would need to work with en- 
semble averages over a large number of realizations, as well as to employ spectra in 



300 


M. T. Landahl, J. Kim and P. R. Spalart 



0 


.5 


1.0 


y/5 



FIGURE 2. Profiles of the root-mean-square values of q in Eqn. (1). 


a frame of reference converted with the local mean velocity. If found to be sound, 
this model will lead to a reasonably simple procedure for determining the Reynolds 
stresses and other statistical quantities through a comparatively simple linear calcu- 
lation making use of a universal model for the nonlinear processes in the near-wall 
region, the statistics of which may be found from numerical simulations carried out 
at modest Reynolds numbers. 
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